Ultimate-state scaling in a shell model for homogeneous turbulent convection 
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An interesting question in turbulent convection is how the heat transport depends on the strength 
of thermal forcing in the limit of very large thermal forcing. Kraichnan predicted [Phys. Fluids 5, 
1374 (1962)] that the heat transport measured by the Nusselt number (Nu) would depend on the 
strength of thermal forcing measured by the Rayleigh number (Ra) as Nu ~ Ra^^^ with possible 
logarithmic corrections at very high Ra. This scaling behavior is taken as a signature of the so- 
called ultimate state of turbulent convection. The ultimate state was interpreted in the Grossmann- 
Lohse (GL) theory [J. Fluid Mech. 407, 27 (2000)] as a bulk-dominated state in which both the 
kinetic and thermal dissipation are dominated by contributions from the bulk of the flow with the 
boundary layers either broken down or playing no role in the heat transport. In this paper, we study 
the dependence of Nu and the Reynolds number (Re) measuring the root-mean-squared velocity 
fluctuations on Ra and the Prandtl number (Pr) using a shell model for homogeneous turbulent 
convection where buoyancy is acting directly on most of the scales. We find that Nu~ Ra^'^^Pr^'^^ 
and Re~ Ra^'^^Pr"^''^, which resemble the ultimate-state scaling behavior for fluids with moderate 
Pr, but the presence of a drag acting on the large scales is crucial in giving rise to such scaling. This 
suggests that if buoyancy acts on most of the scales in the bulk of turbulent convection at very high 
Ra, then the ultimate state cannot be a bulk-dominated state. 

PACS numbers: 47.27.-i, 47.27.te 



I. INTRODUCTION 

In Rayleigh-Benard convection, fluid confined in a box 
is heated from below and cooled on top. When the tem- 
perature difference is large enough, convective motion 
sets in. The flow state is characterized by the geome- 
try of the box and two control parameters: the Rayleigh 
number (Ra), which measures the strength of the thermal 
forcing and the Prandtl number (Pr) , which is the ratio of 
the diffusivities of momentum and heat of the fluid. The 
two parameters are defined by Ra = agAL'^ / {vk), and 
Pr = where A is the temperature difference, L is the 
height of the box, g the acceleration due to gravity, and 
a, v, and k are respectively the volume expansion coef- 
ficient, kinematic viscosity and thermal diffusivity of the 
fiuid. When Ra is sufficiently large, the convective mo- 
tion becomes turbulent. Turbulent Rayleigh-Benard con- 
vection has been a system of great research interest (see, 
for example, [l|, Q for a review). In particular, an in- 
teresting question is how the heat transport of the fiuid, 
measured by the Nusselt number (Nu), which is defined 
as the actual heat transport normalized by that when 
there were only pure conduction, depends on Ra in the 
limit of very high Ra. More than 40 years ago, Kraichnan 
predicted [3] that in this asymptotic limit. 



for 0.15 < Pr < 1. Here Reo — uoL/{2u) is the Reynolds 
number measuring the root-mean-squared horizontal ve- 
locity fiuctuations uq at mid-height L/2. Such a scal- 
ing behavior of Nu ~ Ra^^^ and Reo ~Ra^/^ is taken to 
be a signature of the so-called ultimate state of turbu- 
lent convection. This predicted asymptotic increase of 
Nu as Ra^/^ is stronger than the observed dependence 
of Ra''' with 7 around 0.3 at moderate Ra. According to 
Kraichnan, the convective eddies, produced in the bulk 
of the turbulent convective fiow, generate turbulent shear 
boundary layers near the walls and it is the small-scale 
turbulence present in these boundary layers that domi- 
nates and enhances the heat transport. Thus, the shear 
boundary layers play a crucial role in giving rise to the 
ultimate-state scaling in Kraichnan's work. 

On the other hand, a recent theory proposed by Gross- 
mann and Lohse [3] agrued that the kinetic and thermal 
boundary layers would either break down or do not con- 
tribute to the energy and thermal dissipation and thus 
do not play any role in the heat transport at very high 
Ra. In this bulk-dominated state, the Grossmann-Lohse 
(GL) theory predicted that for fluids with moderate Pr: 



Nu - Rai/2pj.i/2 
Rcisc - Rai/2pr-i/2 



(5) 
(6) 



Nu - Ra^/2(lnRa)-3/2pri/2 
Reo ~ Rai/2(lnRa)-i/2pr-i/2 



for Pr < 0.15 and 



Nu ^ Rai/2(lnRa)-3/2pr-i/4 
Reo ~ Rai/2(lnRa)-i/2pr-3/4 



(1) 
(2) 



(3) 
(4) 



where ^e^sc = UL/u is the Reynolds number measuring 
the mean large-scale circulating fiow velocity U near the 
boundaries. Thus GL predicted the same scaling of Ra^/^ 
for Nu and Rcl^c for fluids with moderate Pr. The Pr- 
dependence predicted in the GL theory agrees with that 
of Kraichnan for Pr < 0.15 but not for 0.15 < Pr < 1. 
We emphasize that although the predicted asymptotic 
dependence of Nu and Rclsc or Reo are both Ra^/^ in 



2 



the two theories, the assumed roles of the boundary layers 
in heat transport in the asymptotic regime are rather 
different. 

The ultimate state of turbulent convection has been 
elusive Q in that definitive experimental evidence is lack- 
ing. An increase in the Nu-Ra scaling exponent was 
found around Ra=10^^ by Chavanne et al. p, 0| in ex- 
periments using low temperature helium gas, and was 
interpreted as the transition to the ultimate state. How- 
ever, similar experiments by Niemela et al. ^] showed 
that the measurements of Nu can be well represented by 
Nu ^ Ra° "^°^ for Ra up to around 10^^. This puzzling 
discrepancy in the two experiments remains unresolved. 
The situation is further complicated by the increasing 
difficulty to keep the experiments within the Boussinesq 
approximation at high Ra @. 

GL's work led to the idea of attaining the ultimate- 
state scaling at moderate Ra by an artificial destruction 
of the boundary layers [Toj . By numerically simulating 
the bulk of turbulent Rayleigh-Benard convective flow, 
modeled by three-dimensional homogeneous turbulent 
convection with periodic boundary conditions flU , Lohse 
and coauthors [13, [13 reported results that are consis- 
tent with Eqs. (O and (O, when Re^sc is replaced by 
the Reynolds number measuring the root-mean-squared 
velocity fluctuations. It has been found that [111, [Si- 
for three-dimensional homogeneous turbulent convection, 
buoyancy is relevant only at the larg est scales. On the 



other hand the Bolgiano length [14|, which is an esti- 
mate of the length scale above which buoyancy forces are 
dominant, increases with Ra at moderate Ra. Thus it is 
unclear whether or not buoyancy remains to be relevant 
only at the largest scales at very high Ra. 

It is therefore interesting to study the scaling of heat 
transport using a model for homogeneous turbulent con- 
vection in which buoyancy is acting directly on most of 
the scales. In this paper, we perform such a study us- 
ing a shell model for homogeneous turbulent convection, 
and buoyancy is acting on most of the scales in some pa- 
rameter range. We investigate the scaling behavior of Nu 
and the Reynolds number (Re) measuring the root-mean- 
squared velocity fluctuations. The paper is organized as 
follows. In Secini we describe the shell model for homo- 
geneous turbulent convection used in the present study, 
define Nu, Re and Ra in the model, and derive two exact 
results. We present and discuss our results of the depen- 
dence of Nu and Re on Ra and Pr in Sec. IIIII In Sec. IIVI 
we understand why our observed result of the average en- 
ergy dissipation e is different from that predicted in GL. 
With the relative simplicity of the shell model, we can 
derive analytically results for the scaling behavior of Nu 
and Re. We will illustrate this in Sec. |Vl and also show 
that the theoretical results agree well with the numeri- 
cal observations. Finally, we summarize and conclude in 
Sec.lVll 



II. THE SHELL MODEL FOR HOMOGENEOUS 
TURBULENT CONVECTION 

Homogeneous turbulent convection, which represents 
the bulk of turbulent Rayleigh-Benard convection, has 
been proposed [ll| as a three-dimensional convective flow 
in a box, with periodic boundary conditions, driven by 
a constant temperature gradient along the vertical direc- 
tion. In Boussinesq approximation [l5|, the equations of 
motion read p^ : 



du 

'at 

80 



u-Vu = -Wp + vV'^u + ag9z (7) 
-u-ve = nV'^e + Pu^ (8) 



with V-u = 0. Here, u is the velocity, p is the pressure di- 
vided by the density, 6 = T— {Tq — (3z) is the deviation of 
temperature T from a linear proflle of constant temper- 
ature gradient of —(3, Tq is the mean temperature, and 
£ is a unit vector in the vertical direction. A dynamical 
shell model for this system has been proposed by Bran- 
denburg [l^. Shell model is constructed in a discretized 
Fourier space with fc„ = fco/i", n = 0, 1, . . . , iV — 1, 
being the wavenumber in the nth shell, and h and ko 
are customarily taken to be 2 and 1 respectively. Shell 
models for homogeneous and isotropic turbulence have 
been studied extensively and proved to be successful in 
reproducing the scaling properties observed in experi- 
ments (irj . In Brandenburg's model, the velocity and 
temperature variables m„ and 9n are real and satisfy the 
evolution equations: 

(9) 

dt 



"vk^Un + agOn 

akn{Un-lOn-l - hUrfin+l) 



+ 6fc„(u„6'„_i - hUn+l9n+i) - nk^On + PUn (10) 



where a, b, a, and b are positive parameters. 

For this shell model, it was found [l^ that when b/a 
is larger than some critical value around 2, the effect 
of buoyancy is greater than the average energy dissipa- 
tion rate for most shells. Moreover, buoyancy directly 
affects the statistics of the system such that the scaling 
behavior of Un and 0„ is given by the Bolgiano-Obukhov 



scaling [ll[Ti (i 



kn^^^, On kn'^'") plus corrections 



-1/5 



rather than Kolmogorov 1941 scaling [20] {i 



kn 



1/3 



On ^ kn ^^'^) plus corrections. In other words, buoyancy 
is directly acting on most of the scales when b/a is suf- 
ficiently large. Thus we shall focus on b/a lar ge i n the 
present work. It was reported in earlier studies [la] that 
the value of b/a controls the direction of energy transfer. 
For large b/a, there is an inverse energy transfer from 
small to large scales. The direction of energy transfer 
can be quantified by the sign of the energy transfer rate. 
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FIG. 1: Fu{k„) for a = 0.01, b = a = b = ag = 1, (3 = 100, 
iy = K = 10"*. 



Multiply Eq. Q by u„ , we get 
dE 

= Fu{kn) - Fu{kn+i) - vkl^v?^ + agUnOn (11) 
where £"„ = is the energy in the nth shell and 

Fu ikri) — kji [aUn-l + hUn)Un-lUn (12) 

is the rate of energy transfer or energy flux from (n— l)th 
to nth shell. As shown in Fig.[Tl {Fu{kn)) is indeed neg- 
ative, confirming that, on average, energy is transferred 
from large n (small scales) to small n (large scales) when 
h/a is large. As a result, a drag acting on the largest 
scale has to be added to dampen the growth of energy at 
large scales so that the system can achieve a statistically 
stationary state [2l|, [l^, [1^ • Thus we modify Eq. ([9]) to 



= akn{ul_-^ - hUnUn+l) + 6fc„(M„U„_i 



(13) 



where /uof^n.o? with / > 0, is a linear drag term acting 
only on the first shell n — 0. Also, Eq. (fTT|) becomes 



dEn 
dt 



Fuikn) - Fu{kn+l) 



vklul 



agur, 



(14) 

Next we need to define Ra, Nu and Re in the shell 
model noting that Pr is given by the usual definition of 
v/k. The definitions of Ra and Re are straightforward: 
we only need to replace L by 1/fco such that 



Ra 



Re 



vko 



(15) 
(16) 



As for Nu, we recall its definition in turbulent Rayleigh- 
Benard convection as: 

_ {u,{T-To)-KdT/dz)A _ {u.9)v 

kA/L kA/L ^ ' 

where {• • ■)a denotes an average over (any) horizontal 
plane of the convection cell and time, and {■ ■ ■)v denotes 
an average over the whole volume of the convection cell 
and time. Thus we define 



Nu = " ' + 1 



(18) 



accordingly with (• • • ) denotes an average over time. 

One can derive two exact results in exact analogy to 
those derived in the case of turbulent Rayleigh-Benard 
convection Multiply Eq. (fTO|) by 0„, we get 

7 Q 

= Fe{kn) - Fg{kn+i) - Kklel + /3u„0„ (19) 

dt 

where Sn = ^^i/2 is the entropy, which is proportional to 
the volume integral of temperature fluctuations in Bousi- 
nessq approximation, in the nth shell and 



(20) 



is the rate of entropy transfer or entropy flux from (n — 
l)th to rtth shell. In the statistically stationary state, 
summing Eq. and Eq. over n and using Eqs. 
and (jlSp give the exact results as 

etotai - '^''k^ (Nu - 1) Ra Pr-^ (21) 
X = i^P^ Nu (22) 

where the total energy dissipation rate etotai is given by 

£total = e + £drag (23) 

Here e is the average energy dissipation rate given by 

e^iyY^kliul) (24) 

n 

and Cdrag is the average rate of energy dissipation due to 
the large-scale drag: 

edrag = /(Uq) (^5) 

The average thermal dissipation rate x is defined as 

XEE«^A;2(0+'«/3' (26) 



in accordance with x = «;((Vr)^)y = K((V6')^)y -I- 
for turbulent Rayleigh-Benard convection. 
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III. RESULTS FOR Nu(Ra, Pr) AND Re(Ra, Pr) 

We numerically integrate Eqs. pU]) and using 
fourth-order Runge-Kutta method with an initial con- 
dition of Un = On = except for a small perturbation 
of 6n in an intermediate value of n. We use a = 0.01, 
b = 5. = b = ag = l, P = 100, / = 0.5, N = 30, and 
vary v and k to study Nu and Re as a function of Ra and 
Pr. Five moderate values of Pr ranging from 0.1 to 2 are 
studied. 

The dependence of Nu on Ra for each Pr are shown in 
Fig. [21 It can be seen that for each Pr, Nu scales with 
Ra: Nu = A(Pr)Ra^^ with 71 = 0.500 ± 0.001. The Pr- 
dependence of Nu is found to be: A(Pt) = CiPr'''^ with 
72 — 0.51 ± 0.01, as shown in Fig. [31 Thus we have 




FIG. 2: Dependence of Nu on Ra for for Pr = 0.1 (cir- 
cles), Pr=0.25 (squares), Pr=0.5 (triangles), Pr=l (stars), 
and Pr=2 (diamonds). The solid lines are the linear least- 
square fits of log Nu versus log Ra. 




FIG. 3: The prefactor ^(Pr) as a function of Pr. The solid 
line is the linear least-square fit of log A versus log Pr. 



Similarly, we study the dependence of Re on Ra for 
each Pr. As seen in Figs. H and [5l Re = B(Pr)Ra''3 
and B{Pt) = C2 Pr'^% where 73 = 0.500 ± 0.001 and 
74 = -0.50 ±0.01. Thus 

Re = C2Ra°-^"°±°-°°iPr-°-5°±°°i (28) 

Hence our results for Nu and Re are consistent with 
Eqs. (O and ([B]), when Re^sc is replaced by Re, and 
are also consistent with the numerical results [10, ob- 
tained in the three-dimensional homogeneous turbulent 
thermal convection in which buoyancy only acts on the 
largest scales. 



Pi 




Ra 



FIG. 4: Dependence of Re on Ra for different values of Pr, 
with the same symbols as in Fig. [2 The sofid fines are the 
linear least-square fits of log Re versus log Ra. 



lOi 




0.1 1 
Pr 

FIG. 5: The prefactor i3(Pr) as a function of Pr. The solid 
line is the linear least-square fit of log B versus log Pr. 

Moreover, we note that the Pr-dependence of Nu and 
Re is consistent with GL and not with Kraichnan's re- 
sult for 0.15 < Pr < I at very high Ra. The signifi- 
cance of Nu - (Ra Pr)i/2 and Re ~ (Ra/Pr)!/^ is that 



the heat transport and the root-mean-squared velocity 
fluctuations are independent of v and k. Thus for flu- 
ids with moderate Pr, the scaling Nu ^ (RaPr)^/-^ at 
very high Ra is in line of a 'central dogma in turbu- 
lence' j5] that the effects of turbulence become indepen- 
dent of viscosity and thermal diffusivity when Re is suf- 
ficiently large. On the other hand, the rigorous upper 
bound of Nu< 0.167 Ra^/^ — 1 [24], for convection in a 
layer of fluid with no sidewalls, indicates [l^ that the 
dependence of Nu ^ (Ra Pr)^/^ cannot hold for fluids 
with very large Pr. 

In the derivation of Eqs. ^ and (O in the GL the- 
ory, there are two key intermediate results, which are 
the estimates of e and x in the bulk-dominated regime 
for moderate Pr. In the shell model, L ~ 1/fco, Re^sc 
becomes Re, and these results translate to 



AGL) 



ly^k^ Re^ 
RePr 



Next we investigate the validity of Eqs. 
Fig. [SI we see that 



(29) 
(30) 

and §U\). In 




FIG. 7: x/ {i^0^) ^ a function of RePr with the same symbols 
as in Fig. [J] The solid line is the linear least-square fit to all 
the data points in the log- log plot. 



= C3Re2-48±o.o2 



(31) 



Thus the dependence on Re is approximately Re^/^ in- 
stead of Re"^. On the other hand, as can be seen in Fig. [T] 



J, = C.(RePr) 



1.000±0.001 



(32) 



in good agreement with Eq. ([301 




FIG. 6: e/ (v^k"^) as a function of Re for different values of Pr 
with the same symbols as in Fig. [2] The solid line is the linear 
least-square fit to all the data points in the log-log plot. 



The observed results of Eqs. ([?T|) and ([5^ together 
with the exact results Eqs. (|2ip and ([^^ imply that ttotai 
cannot be dominated by e otherwise we would have Re 




FIG. 8: edrag / iy'^ ko) as a function of Re with the same sym- 
bols as in Fig. [2] The solid line is the linear least-square fit 
to all the datapoints in the log-log plot. 



(Ra/Pr)i/('^"i) K, (Ra/Pr)2/3, which is in contradiction 
to our observed dependence of Re ~ (Ra/Pr)^/^. Thus, 
^totai has to be dominated by edrag- Moreover, tdrag has 
to scale as Re'^. Indeed we find that e ^ etotai such that 



^total ~ ^drag 

and as shown in Fig. [51 



CsRe' 



3.00±0.01 



(33) 



(34) 



giving consistent results as expected. 

Thus the observed scaling results of Nu and Re 
[Eqs. ([77)1 and ([^5|l ] depend crucially on the presence of 
a large-scale drag. Such a damping mechanism at the 
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largest scales cannot exist by itself in the bulk of turbu- 
lentthermal convection but could be resulted from the in- 
teraction of the boundaries with the buoyancy-generated 
inverse transfer of energy from small to large scales. This 
suggests that when buoyancy is acting directly on most 
scales, the ultimate state, if exists, cannot simply be a 
bulk-dominated flow state. Instead the boundaries must 
play a crucial role. 



IV. UNDERSTANDING THE DEPENDENCE 
OF e ON Re 

In this section, we discuss how one can understand the 
observed dependence of e on Re approximately as Re^^^. 
The average energy dissipation rate in each shell increases 
with n up to a maximum at the dissipative scale, whose 
shell is denoted by the shell number rid, then decreases 
again. Thus e can be approximated as : 



(35) 



where Di is a number of order 1. The dissipative 



wavenumber /c„^ can be estimated as usual as 



(36) 



with D2 being a number of order 1. Now (u^J has good 
scaling behavior in fc„ [H, for < rt < ti^: 



-2>) 



Moreover, 



(37) 



(38) 



where D3 is a number of order 1. Putting Eqs. 
together, we get 



,31,4 



Rei+'; 



(39) 



As discussed, when buoyancy is acting on most of 
the scales, the scaling behavior is given by Bolgiano- 
Obukhov plus corrrections, thus 77 ~ 3/5. This gives 
4/(1 -I- 7y) « 5/2, as observed [see Eq. ^H]. On the other 
hand when buoyancy is only acting on the largest scales, 
Un obeys Kolmogorov 1941 scaling plus corrections p^ . 
In this case, 77 « 1/3, which leads to 4/(1 + 77) « 3, giving 
the same intermediate result used in the GL theory [see 
Eq. (I29p ]. Such a dependence of Re^ can be easily un- 
derstood as follows. When buoyancy is only acting as 
a driving force at the largest scales, we have the usual 
energy cascade. From Eq. (jlip . we have 



using Eq. Using Eq. Eq. ^ leads immedi- 

ately to e ^ ly^kQ Re"^. When buoyancy is acting on most 
of the scales, there is, however, no longer a cascade of 
energy but only a cascade of entropy 18]. In particular 
(F„(fc„^)) ^ (F„(fci)) and, as a result, e docs not scale 
as Re"^. In other words, the result Eq. (|29|) estimated in 
the GL theory does not hold when buoyancy is directly 
acting on most of the scales with the dynamics governed 
by an entropy cascade. 



V. ESTIMATES OF x AND earag AND THUS 

Nu(Ra, Pr) AND Re(Ra, Pr) 

In the shell model, there is always an entropy cascade 
with {Fg(kn)) ~ X being independent of fc„ in the inter- 
mediate scales. Thus 

X^{Feik^))^ko(uiy/'{el) (41) 

using Eq. (^0]) . On the other hand, Eq. implies that 

Fe{k,) « (3{uo0o) « /3(u2)V2(02^i/2 (42) 

Comparing Eqs. ([4T|) and ([42]). we get 

ko{diy/^^l3 (43) 



Putting Eq. gSl) into Eq. (gTI) and using Eq. ([38]), we 
obtain 



X •■ 



0- 



/2 



(44) 



(45) 



which leads immediately to 

X « K0^Re Pr 

as observed [see Eq. (15^ ]. 

From Eqs. and ([351), ^e have 



edrag ^ fklv^^e" = v''4^L= '^Re^ (46) 

For €drag/ (v^kQ) to depend only on Ra and Pr, we require 



/ = /i(Ra, Pr)Vag^ 



(47) 



for some function h of Ra and Pr when Ra and Pr are 
varied. In our numerical calculations, ag, P and / are 
kept fixed while v and n are varied to get different Ra 
and Pr. This amounts to taking 



/i(Ra,Pr) = ho 



(48) 



for some fixed constant Hq. Using Eqs. (PTjl . ([551) . iH^ - 
([47|. and we get 



(40) 



Nu Ra^/^ w hoP^^^Re^ 



(49) 
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On the other hand, Eqs. and ([45]) imply 

Nu « Re Pr (50) 

Solving Eqs. dH]) and we get 

Nu « -LRai/2pri/2 (g^) 

Re « J_Rai/2pi.-i/2 (52) 
ha 

which are exactly the ultimate-state scaling observed [see 
Eqs. (EZl) and Substituting Eqs. (gT]), (gS]), and 

([5^ into Eq. one gets 

edrag « i^'A^^/lgRe^ (53) 

which is in good agreement with our observation [see 
Eq. jSl]. 

We have demonstrated that the presence of a drag that 
acts on the largest scales is crucial for the observation of 
the ultimate-state scaling of Nu and Re. In our calcula- 
tions and derivation, we employ a linear drag of the form 
fuQ. An immediate question that arises is whether the 
scaling laws of Nu and Re depend on the specific mathe- 
matical form of the large-scale drag used. To answer this 
question, we replace fuoSn,o in Eq. ((T3|) by the general 
form of a nonlinear drag: fu^~^Sn,o, where m > 2 is an 
integer, and repeat our analysis to get a new estimate 
of the generalized average energy dissipation rate due to 
the drag, idrag- Now 

edras = no « f{ulr'^ (54) 
Using Eq. dSH]), we get 

edrag ~ fc^^ /iQ (RaPr" ^ ) ^ Re™ (55) 

where 

/ = hoiy/^)'-"'k^-^ (56) 

for some fixed value Hq as Ra and Pr are varied. Using 
^totai ~ edrag, Eqs. ([2l]) and ([55]), we now have 

Nu Ra(™-i)/2 w /loPr'^+i'/'Re™ (57) 

in place of Eq. ([19]). Together with Eq. ^0^, which re- 
mains the same, we get 

Nu « /io-'/^"-')Rai/2pri/2 (58) 
Re « /io"'/'""'^Ra'/'Pr^'/' (59) 

In other words, the scaling results of Nu ~ (RaPr)^/^ and 
Re ~ (Ra/Pr)^/2 remain valid for the general large-scale 
drag of /M™~^(Sn,o, for to > 2. 



VI. CONCLUSIONS 

An interesting question in turbulent thermal convec- 
tion is how Nu and Re depend on Ra and Pr at very high 
Ra. Both the theories by Kraichnan Q and Grossmann 
and Lohse [3| predicted that in this limit, Nu and Re 
would scale with Ra^/^ foj- fluids with moderate Pr. This 
kind of scaling behavior is taken to be the signature of the 
ultimate state of turbulent convection. However, the two 
theories have rather different assumptions about the role 
of the boundary layers in heat transport. According to 
Kraichnan, convective eddies produced in the bulk gener- 
ate turbulent shear boundary layers and the turbulence of 
which enhances heat transport. Thus we interpreted that 
in this picture of Kraichnan, buoyancy is directly acting 
on all scales in the bulk and that the boundary layers play 
an important role in heat transport. On the other hand, 
Grossmann and Lohse argued that in the limit of high 
Ra, the boundary layers would either break down or not 
contribute to the energy and dissipation and thus play 
no role in heat transport. Studying numerically three- 
dimensional homogeneous turbulent convection in which 
buoyancy is ac ting only on the largest scales, Lohse and 
coauthors [l^, Il2| reported scaling behavior of Nu and 
Re that is consistent with the ultimate-state scaling. 

In the present work, we have studied the scaling behav- 
ior of Nu and Re using a shell model of homogeneous tur- 
bulent convection in which buoyancy acts on most of the 
scales. In this model, buoyancy modifies the statistics of 
the velocity fluctuations such that the statistics are given 
by Bolgiano-Obukhov plus corrections instead of Kol- 
mogorov 1941 plus corrections [1^1 . Moreover, there is an 
average inverse energy transfer from small to large scales 
such that a large-scale drag has to be present for the 
system to achieve statistical stationarity. Such a large- 
scale drag cannot exist by itself in the bulk of turbulent 
convection but could be resulted from the interaction of 
the inverse energy flow with the boundaries. Thus, when 
buoyancy is acting directly on most of the scales in the 
bulk of turbulent convection, the boundary layers would 
play a crucial role and, as a result, the flow cannot be 
bulk-dominated. In this case, we have found that the de- 
pendence of Nu and Re on Ra and Pr is again consistent 
with the ultimate-state scaling. 

Because of the relative simplicity of the shell model, 
we can understand analytically the scaling behavior of 
Nu and Re. The two exact results Eqs. and ((^ 
derived for the shell model are in exact analogy to those 
derived for turbulent Rayleigh-Benard convection. As is 
clearly demonstrated by the GL theory, to get the Nu 
and Re scaling, the whole task is to estimate Ctotai and 
X- In the shell model, there is always a cascade of en- 
tropy so that X is given by k/S^ RePr [see Eq. (g5)l ] . When 
buoyancy is acting on most of the scales, etotai is dom- 
inated by €drag- For a linear drag, edrag is estimated 
as ly^k^ (Ra/Pr)i/2Rg2 jggg ^^^^ ^^lat in 

this case, e is given by j/^fc^Re^/^ [Eq. ^ with r] w 3/5 
when buoyancy acts on most of the scales] rather than the 



8 



prediction of i/^k^Ke'^ by the GL theory. Putting these 
estimates into the two exact results, we get the uhimate- 
state scahng of Nu - (RaPr)!/^ and Re ~ (Ra/Pr)!/^. 
On the other hand, when buoyancy is acting only as a 
driving force on the largest scales, edrag is negligible and 
etotai given by e as usual. In this case, the statistics of the 
temperature resemble those of a passive scalar, and e is 
given by v^^k^Re^ [Eq. ([55]) with 77 « 1/3 when buoyancy 
acts only on the largest scales] , the usual result obtained 
in inertia-driven turbulence without buoyancy, which is 
also the result derived in the GL theory. In this case 
Eqs. and ([50]) hold, leading to the ultimate-state 
scaling as shown in the GL theory. 

Hence there are two different physical scenarios 
that can give rise to the ultimate-state scaling of 
Nu - (RaPr)i/2 and Re ~ (Ra/Pr)!/^ for fluids with 
moderate Pr. In the first scenario, which is illustrated 
in the present work, buoyancy is acting directly on most 
of the scales of the bulk of turbulent convection and, on 
average, energy is transferred from small to large scales. 
An effective damping at the largest scales, which can be 
provided by the interaction of the inverse energy trans- 
fer with the boundaries, is crucial. In the second sce- 



nario, buoyancy is acting only as a driving force on the 
largest scales, temperature in the bulk of the convective 
flow is behaving like a passive scalar statistically, and the 
boundary layers play no role in heat transport. The first 
scenario is in accord with the physical picture presented 
in Kraichnan's work Q while the second scenario is in ac- 
cord with that proposed by the GL theory f3|. The next 
question would be whether or not buoyancy acts directly 
on most of the scales in the bulk of turbulent Rayleigh- 
Benard convection at very high Ra, and the answer of 
which would help to distinguish which scenario is physi- 
cally relevant. 
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